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Abstract 

Many new applications of ferroelastic materials require a better understanding 
of their dynamics that often involve phase transformations. In such cases, an impor- 
tant prerequisite is the understanding of wave propagation caused by pulse-like load- 
ings. In the present study, a mathematical model is developed to analyze the wave 
propagation process in shape memory alloy rods. The first order martensite trans- 
formations and associated thermo-mechanical coupling effects are accounted for by 
employing the modified Ginzburg-Landau-Devonshire theory. The Landau-type free 
energy function is employed to characterize different phases, while a Ginzburg term 
is introduced to account for energy contributions from phase boundaries. The effect 
of internal friction is represented by a Rayleigh dissipation term. The resulted non- 
linear system of PDEs is reduced to a differential-algebraic system, and Chebyshev's 
collocation method is employed together with the backward differentiation method. 
A series of numerical experiments are performed. Wave propagations caused by im- 
pact loadings are analyzed for different initial temperatures. It is demonstrated that 
coupled waves will be induced in the material. Such waves will be dissipated and 
dispersed during the propagation process, and phase transformations in the material 
will complicate their propagation patterns. Finally, the influence of internal friction 
and capillary effects on the process of wave propagation is analyzed numerically. 
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1 INTRODUCTION 

In the past decades, different aspects of Shape Memory Alloys (SMA) have been inves- 
tigated intensively by mathematicians, physicists, and engineers [S]. This interest in a 
larger scientific community is due to SMA unique properties of being able to convert 
thermal energy into mechanical and vice versa. These properties are promising for many 
applications of SMAs, including mechanical and control engineering, biomedicine, com- 
munication, robotics to name just a few [B]. Motivated by application developments of 
advance composite materials involving SMAs, nonlinear wave propagations in these ma- 
terials have been investigated as a stepping ground for the prediction and understanding 
of dynamic response of the composite under various dynamic loadings [H SI fT3] . 

Compared to wave propagations in conventional solid materials, the impact induced 
wave propagations in materials such as SMAs requires delicate treatments as additional 
difficulties arise due to phase transformations [U HI [I3]. In general, impact loadings of 
these materials will cause nonlinear thermo-mechanical waves which are similar to those of 
other thermo-elastic materials under impact loadings. The difference of wave propagation 
in conventional solids and those in the ferroelastic materials such as SMAs is that the first 
order martensitic transformation may be induced in the latter case. The transformation 
is reversible, and its native nonlinearity and hysteresis will have a substantial influence 
on the wave propagation and will make the wave propagation patterns more complicated 

The first step to the modelling of impact induced wave propagations and phase trans- 
formations is a sound constitutive theory upon which the entire model can be built p!l[T3]. 
Various constitutive models have been proposed on mesoscale and microscale to capture 
the phase boundary movement induced by the dynamical loadings [21 12] • Example of 
such a constitutive model can be found in Ref. [T[ |2] where an one- dimensional model 
for the modelling of shock wave propagations with phase transformation was constructed 
on the basis of a non-convex Helmholtz free energy function. Under this approach the 
entire structure was split into different domains due to the phase transformation and 
the movement of boundaries between the domains was modelled using the "jump condi- 
tions". This approach is suitable for microscopic problems, while for many engineering 
applications a model is required at macroscale. In Ref. [101 [E], the dynamic behaviour 
of phase boundaries was modelled using a thermo-mechanical coupling approach. The 
model was based on a linearized constitutive theory, and hence its application potential 
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was inherently limited. 

For many engineering applications, the dynamic response of SMA materials caused 
by impact loadings need to be better understood at macroscale for design and control of 
SMA-based devices. For this purpose, displacement and temperature evolutions in the 
material are normally sought. Models at mesoscale may not be sufficient for this purpose 
as another model needs to be constructed to link macroscale properties and mesoscale 
domain structures. Another aspect of modelling the dynamics of such materials as SMAs 
under impact loadings is the thermo-mechanical coupling effects. In most of the existing 
investigations, the thermal dynamics are either neglected [Hill [12], or modelled separately 
from the mechanical dynamics [IHIIIZ]. However, the physics-based models should account 
for the intrinsic coupling of thermal and mechanical fields in SMAs. When the SMAs are 
used for damping purposes or for other purposes where the conversion of energy between 
the thermal and mechanical fields is essential, the coupling effects are expected to be 
particularly important, and the constitutive theory should be constructed by taking into 
account both fields simultaneously. 

In this paper, the nonlinear thermo-mechanical wave propagations in SMA rods in- 
duced by impact loadings are modelled and analyzed at macroscale. To capture the 
thermo-mechanical coupling and nonlinear nature of the phase transformations, the Ginzburg- 
Landau- Devonshire theory is applied for the modelling of the nonlinear dynamics. The 
governing equations for the mechanical field are obtained by minimizing the mechani- 
cal energy, while those for the thermal field are obtained by using the conservation law 
of internal energy. The intrinsic coupling of the two fields is built-in into the model by 
including both fields in the potential energy functional. In the following sections, a mathe- 
matical model describing SMA dynamics is developed based on a system of coupled partial 
differential equations which is re-cast in the form of differential algebraic equations, and 
Chebyshev' collocation method is employed together with the backward differentiation 
formula to integrate the resulting system. Nonlinear wave propagation patterns caused 
by an impact stress loading at one end of the rod are simulated with different initial tem- 
peratures and computational parameters. Finally, the influence of phase transformation 
on the wave propagations is analyzed numerically, along with the influence of other effects 
such as internal friction and capillary effects. 



2 THE INITIAL-BOUNDARY VALUE PROBLEM 

We restrict our analysis to one-dimensional cases as sketched in Figure (jlj). The SMA 
rod under consideration occupies an interval [0, L], and is subjected to an impact loading 
from the right end x = L, while the other end x = is fixed. The rod is thermally 
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insulated at both ends so there is no heat loss (gain) to (from) the ambient environment. 
Under external loadings, a material point x in the SMA rod will be moved to a new 
position X + u{x,t) due to deformation, where u{x,t) is the longitudinal displacement at 
time t. Function u{x,t) is assumed to be continuous at any time t and position x based 
on the continuity of the rod at macroscale. The stress a is related to the deformation e 
by a{x,t) = Af{e{x,t)) where e{x,t) = du{x,t)/dx is the strain. 

By now, it is well understood that the first order phase transformation in SMAs holds 
the key to unique properties of the material, such as the shape memory and pseudo- 
elastic effects. Therefore, it is expected that the adequate mathematical model for the 
dynamics of SMAs at macroscale should be able of capturing the first order martensite 
phase transformations. On the other hand, the intrinsic coupling of mechanical and 
thermal fields should also be captured by the model. To satisfy these requirements, the 
mathematical model based on the modified Ginzburg-Landau-Devonshire theory has been 
established [3 El [13]: 

where fci, /c2, fcs and kg are material-specific constants, v internal friction coefficient, p the 
density of the material, 9i the reference temperature, the specific heat capacitance, k 
heat conductance, and / and g mechanical and thermal loadings, respectively. 

It is essential that the above model is constructed on the basis of the potential energy 
function JF, which a non-convex function of the chosen order parameters and temperature 
9, at mesoscale, according to Ginzburg-Landau-Devonshire theory. It is a sum of local 
energy function [J^i) and non-local energy function {Tg). For the current one-dimensional 
problem, the strain e{x,t) is chosen as the order parameter, and the local free energy 
density can be constructed as the Landau free energy density TiiO, e), while the non-local 
part can be formulated as a Ginzburg term TgiVe): [31 [9l 113] : 

^(^,5) = ri{e,e)+rg{Ve), 

HO.s) ^ M£^,. + |,. + ^,a^ (3) 
Bf 

Tg{Ve) = kg{-)\ 

The local minima of the local term (Landau free energy function) are introduced to 
characterize martensite variants, while the non-local term (Ginzburg term) above accounts 
for inhomogeneous strain field, which represents energy contributions from domain walls 
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and other boundaries among different phases. It will be translated into capillary effects 
at macroscale. 

In order to account for the internal friction, accompanying wave propagations and 
phase transformations, which will be translated into viscous effects at macroscale, a 
Rayleigh dissipation term is also included in the above model by using the following 
term [2]: 

- (4) 

The constitutive relations for the material in the above model at macroscale can be 
obtained by using the thermodynamic equilibrium conditions: 

where T-C{6,e) = T — c^QXwQ is the Helmholtz free energy function. The mechanical and 
thermal fields are intrinsically coupled since the internal energy e is associated with the 
same potential energy as above, and the governing equations for the thermal field are 
formulated using the conservation law for the internal energy [TJ [13]. 

The difficulties associated with the above model are understood better by analyzing 
the profiles of the non-convex potential energy F;, its temperature dependence, and the 
non-convex constitutive curves, as sketched in Figure ([T]). At low temperature Q = 201 
K, as sketched in the top row, there are two local minima in the potential function, which 
are introduced to characterize martensite plus and minus, respectively (in ID cases). The 
stress-strain relation of the material under dynamical loadings will not follow the sketched 
constitutive relations exactly. Instead, there will be jumps from point {A to B) or vice 
versa. Such jumps are associated with the transition from martensite plus to minus or 
vice versa. This is the origin of mechanical hysteresis which will dissipate mechanical 
energy quickly by converting it into thermal form due to the thermo-mechanical coupling. 
The amount of energy converted into thermal form can be estimated by using the area 
enclosed under the hysteresis loop, as marked by the dashed lines in the stress-strain plot 
on the left. When the SMA rod temperature is intermediate {6 = 240 K), there are still 
hysteresis loops, and the jump phenomena become more complicated, since there are three 
local minima present in the potential energy function, which means that austenite and two 
martensite variants may co-exist in the material. The dissipation of mechanical energy 
will be slower in this case since the hysteresis loop become smaller. At high temperature 
{9 = 310 K), it is shown that there are no jump phenomena any longer, therefore no 
hysteresis, because there is only one local minimum. In this case, there will be only 
austenite present and the dynamics become fairly simple. 
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In the current paper, the mechanical loading is implemented in terms of impact stress 
at one of the rod ends. Hence, it is convenient to keep the constitutive relation as an 
extra equation for the model and consider the stress as a dependent variable. This rep- 
resentation will make the treatment of boundary conditions much easier for the current 
discussion. The resulting system of Differential Algebraic Equations (DAE) can be written 
as follows [TS]: 

du dv da d dv , d d'^e 

'dt 

(6) 



^„dt dx~^ dx dx ^ dx dx"^ ' 
de d^e , , de 

''"m=^dx^^^'^'W 
a = ki{e- 61)8 + k2e^ + fcge^ 



where v is the velocity. The mechanical and thermal loadings, / and g, are all set to zero, 
so that only boundary loadings will be taken into account in the current investigation. 

In order to investigate the thermo-mechanical wave propagations in the SMA rod, 
the following boundary conditions are employed for the mechanical and thermal fields 
similarly as in Ref . [3, [2^ : 

de 

dx ^ 
a(L,t) 
d'^u 



de 

dx 



0, 



x=0 
M(0,t) = 0, 

d'^u 



0, 



5(t), 



(7) 



dx"^^ 



0, 



dx"^ 



0. 



where S{t) is a given function describing the stress impact profile. 



3 Wave Propagations 

Since mechanical responses caused by external loadings in most materials are normally 
much faster than thermal ones, which are also more interesting for numerical investigations 
and many applications, the emphasis of the current discussion is put on mechanical waves 
caused by mechanical loadings. 

3.1 Temperature Dependence 

For the analysis of elastic waves in the SMA rod. Equation ([1]) is firstly linearized at the 
point (ctl. El) where a^ = 0: 

(9 , d de ^ d^u 

where the external mechanical loading is dropped out. is the stiffness constant for the 
linearized system, which is temperature dependent. 
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When the SMA is at high temperature {6 = 310 K), only austenite is stable and there 
is no phase transformation, it can be easily calculated that = so can be simply 
formulated as: ki = ki{6 — 6i). 

At lower temperature {6 = 210 K), the stress strain relation is not a monotone curve 
any longer, and the linearization has to be carried out by using at least three points, as 
indicated by the plot in the top row of Figure ([T]) (the three intersections between the 
horizontal axis with the a — e curve). The central one is = 0, which is not a stable 
equilibrium state of the system, and is not interesting for the analysis. The other two 
intersections can be easily calculated using the following condition: 

k^{9 ~ 9i) + k2e^ + kse'' = 0, (10) 

which gives the following formulation (using parameter values given in section 5): 

^,_-W5ZWEM. ,, = ±0.115. (11) 

These two values are associated with strains for martensite plus {el = 0.115) and minus 
{el = —0.115), respectively. The linearized coefficient then can be calculated as follow: 

kL = ki{e-ei) + k2El + ksEl (12) 

which gives ki the same value at the two points with different el values, due to the 
symmetry property. The above analysis indicates that the linearized wave motion in the 
material at martensite plus state will be the same as those at martensite minus state. 

For the cases where the SMA rod temperature is intermediate, the dependence of ki 
on temperature is more complicated since both austenite and martensite might occur 
and there might be 4 values for El to satisfy cr^ = 0, as indicated in the middle row in 
Figure ([1]). Although the symmetry property still exists for martensitic variants, wave 
motions are different between austenite and martensite states. 

Let us consider the following wave solution to the linearized wave motion given by 
Equation ©: 

u = u{x — Vt) = u{z), z = X — Vt, (13) 

where V \s a. constant stands for wave velocity. By substitution, the following relation 
can be easily obtained: 

^^S-l^.-Pn-^t'l. (14) 

where Kl — pV"^ can be positive or negative. The problem is now formulated as an 
ordinary differential equation and its general solution can be written as: 



,(z) = (Ci + C2z)eV(^^-^^')A«, (15) 
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where Ci, C2 are coefficients to be determined by boundary conditions. The viscous term 
is temporarily ignored. 

If the wave velocity V is less than the velocity of sound in the material Vs = {ki/ pf', 
then Kl — pV^ > 0, and there is no limited solution exists. It means that no waves 
can propagate in the material with velocity V smaller than the velocity of sound in the 
material. Since the velocity of sound of the material is temperature dependent, so the 
allowed speed for waves to propagate in the SMA rod varies along with the variation of 
its temperature. 

When the viscous term is also taken into account, then the wave propagation will 
always be accompanied by dissipation effects, which can be characterized by the following 
exponential function: 



where the initial amplitude of the considered waves are assume to be 1. The dissipation 
effects can be estimated by the exponential coefficient = For larger V, faster 
dissipation will be induced. The dissipation effects are independent of the material tem- 
perature. 



3.2 Effects of Ginzburg's Term 

In the wave equation, the term kg^ is resulted from the interfacial energy contribution 
to the potential energy function. It is also called Ginzburg's term, and accounts for the 
capillary effects [3, El [131 E] • It is easy to see that this term is related to dispersion of 
wave propagations. For the analysis, the following solution to the linearized wave equation 
is considered: 

u = sm—{x-Vt), (17) 
A 

where A is the wave length. By substitution, the following relationship involves the wave 
speed V can be obtained if the viscous effects are ignored: 

pV^' = A;l + ^, (18) 

which indicates that the wave propagation speeds are increased due to the contribution of 
non-local potential energy (capillary effects). The dispersion effects caused by non-local 
contributions are stronger for those waves with smaller wave lengthes. If the SMA rod is 
at lower temperature, the linearized stiffness constant kj^ will be smaller, which will make 
the last term in Equation (1181) more pronounced. 
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4 NUMERICAL METHODOLOGY 

As mentioned in the previous section, the development of the numerical methodology 
for simulation of the wave propagation based on the above mathematical model is not a 
trivial task. In particular, given that both dispersion and dissipation of wave propagations 
are present in the physics of the problem, the numerical algorithm for the problem has 
to be able to take care of both dissipation and dispersion numerically, and the accuracy 
of the algorithm will be affected by their treatment. In the present paper, a multi- 
domain decomposition method combined with the Chebyshev collocation methodology is 
the method of choice in addressing the above issues. The compromise made here is that a 
spectral method is employed to take advantage of its better convergence property, while 
the domain decomposition method is chosen for the purpose to reduce the order of basis 
functions for the spectral method when the total node number is large. 

4.1 Chebyshev's Collocation Method 

For the Chebyshev pseudo-spectral approximation, a set of Chebyshev points {xj} are 
chosen along the length direction as follows: 

l-cos(-)j/2, z = l,2,...,Ar. (19) 

Using these nodes, u, v, 9, and a distributions in the rod can be expressed in terms of 
the following linear approximation: 

N 

f{x) = Y,hU^), (20) 
1=1 

where f{x) stands for any of u, v, 9, or a, and fi is the function value at Xj. (f)i{x) is the 
i*'* interpolating polynomial which has the following property: 



U^^) = J . /.' (21) 




It is easy to see that the well-known Lagrange interpolants satisfy the interpolating 
requirements. Having obtained f{x) approximately, the derivative df{x)/dx can be easily 
obtained by taking the derivative of the basis functions 4>i{x) with respect to x: 

dx~^^' dx ' ^^^^ 
1=1 

and similarly for the higher order derivatives. All these approximations can be formulated 
in the matrix form, for the convenience of programming. 



9 



4.2 Multi-Domain Decomposition 

It is known that the spectral methods are able to give a higher accuracy with the same 
number of discretization nodes, compared to finite difference methods or finite element 
methods. On the other hand, when the solution to the problem does not have higher- 
order derivatives, the spectral methods may lead to artificial oscillations due to the Gibbs 
phenomenon. This may be expected for the current problem when the impact induced 
wave propagation is analyzed numerically. To avoid this, a multi-domain decomposition 
method is employed. 

The entire computational domain V = [0, L] is evenly decomposed into P inter- 
vals (subdomains) , with an overlap region between each pair of consecutive intervals, 
as sketched in Figure ([3]): 

p=p 
p=i 



V 



(23) 



where the number of subdomains P is chosen according to the specific problem under con- 
sideration. In each interval, the Chebyshev collocation method discussed in the previous 
section is employed to approximate the solution and its derivatives. 

The coupling between each pair of consecutive intervals can be implemented by setting 
the following requirements: 



p+1' 



(24) 



where the subscript p stands for the interval number, while the superscript n stands for 
the node number in each interval. Variable is the function value at point x'^ (the nth 
node in the pth interval), which could be any of the dependent variables we are solving 
for. Point x^ is actually the same node as Xp_^_^, and x^"^ is the same node as x^^^. 

The derivatives of functions in the overlapped nodes are approximated by taking the 
average of their values evaluated from the two intervals involved: 



dx 



dy_ 

dx 




i=0 

N 



dx 
d(f)i{x) 



(25) 



dx 



>+i/ 



r2 



The approximation to the second order derivatives can be found using the same average 
for the nodes in the overlapped region. 
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4.3 Backward Differentiation Formula Method 



By employing the multi-domain decomposition methods combined with the Chebyshev 
collocation methodology, the given set of partial differential equations Eq. ([6]) can be 
converted into a DAE system, which can be generically written in the following form: 

M— + N{t,X,git)) = 0, (26) 

where X is the vector collecting all the variables we are solving for, M is a singular matrix, 
A?" is a vector collecting nonlinear functions produced by spatial discretization. The 
resultant DAE system is a stiff system and has to be solved by an implicit algorithm. Here 
the second order backward differentiation formula method is employed for this purpose. 
By discretizing the time derivative using the second order backward approximation, the 
DAE system can be converted into an algebraic system at each time level, which can 
formally written as follows: 

M Qx" - 2X"-i + + AtN {tn, X\ g{K)) = 0, (27) 

where n denotes the current computational time layer. For each computational time layer, 
iterations must be carried out using Newton's method for X" by use of X""-*- and X"^^. 
Starting from the initial value, the vector of unknowns X can be solved for at all specified 
time instances employing this algorithm. 



5 Numerical Experiments 

A series of numerical experiments have been carried out to investigate the nonlinear wave 
propagations in the SMA rod involving phase transformations. All experiments reported 
here have been performed on a Au23Cu3oZn47 rod, with a length of 1 cm. The physical 
parameters, except u and kg, for this specific material are taken the same as those in [19] , 
which are listed as follows for convenience: 

ki = 480 g/ms^cmK, k2 = Q x 10^ g/ms'^cmK, fcs = 4.5 x 10^ g/ms'^cmK, 
01 = 208K, p = ll.lg/cm^, = 3. 127 Ag/ms'^cmK, k = 1.9 x 10'"^ cmg/ms^K. 

Experiments indicate that the Ginzburg coefficient kg should be relatively small com- 
pared to ki, so it is taken as kg = lOg/ms^ first, by referring to Ref. [7j. The internal 
friction coefficient is not an easily measurable quantity. In what follows, we assume it 
to be a small fraction (2%) of fci, that is approximately IQg / {cm){ms)). The entire rod 
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is divided into 10 sub-intervals, in each interval there are 15 nodes used for spatial dis- 
cretization. All simulations have been carried out for the time span [0,0.1] ms, and the 
time step-size for the integration is chosen as 2.5 x 10~^ms. 

The first numerical experiment for nonlinear wave propagations in the SMA rod is 
performed with a higher temperature 9 = 310 K, for which there is no phase transfor- 
mation. Other initial conditions are chosen u = v = s = 0, and mechanical loading is 
employed in terms of a stress impact at right end, as follows: 



which can be regarded as an approximation to a pulse stress impact on the SMA rod. 

The mechanical wave propagations are presented by the strain evolution while thermal 
waves by temperature evolution in Figure ([5]). It is shown that the impact induced 
waves start from x = 1 and propagate along the negative x direction, hit the opposite 
boundary at x = and are bounced back. The temperature evolution indicates that there 
are associated thermal waves induced by the stress impact loading, due to the thermo- 
mechanical coupling effects. The propagation patterns of the thermal waves are similar 
to those of mechanical waves. The evolution of the displacement distribution due to the 
stress impact is also presented in Figure ([5]), in the left bottom sub-figure. To clarify the 
patterns of wave propagations, the strain distributions in the SMA rod at three chosen 
time instants are plotted in the right bottom sub-figure in Figure © = 0.01, 0.02, 0.03 
ms, respectively). The arrow attached to each wave profile is to indicate its propagation 
direction. It is seen that the strain distributions in the SMA rod are always smooth and 
no obvious sharp jump occurs, since only austenite is stable with this temperature and 
there is no phase transformation. The amplitude of the wave decreases and the wave 
peak becomes broader during the propagation, it can be easily explained by the fact that 
dissipation effects and dispersion effects, caused by internal friction, capillary effects, and 
thermo-mechanical coupling, are all included in the model. The average wave propagation 
speed can be estimated by the location of the wave frontier or wave peak plotted in the 
figure. With the current initial temperature, the strain wave is bounced back from x = 
and its peak is located around x = 0.65 cm when t = 0.003 ms. This experiment simulates 
the nonlinear thermo-mechanical waves like those in regular thermo-elastic materials. 

The second example deals with the same computational conditions and loading, except 
that the initial temperature now is set 6 = 240K, for which both martensite and austenite 
may co-exist in the SMA rod. The numerical results for this case are presented similarly 
in Figure (EI). It is easy to see that the strain and temperature waves are not as regular 
as those in the first experiment. There are some plateaus clearly shown in the strain 
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and temperature figures, which can be related to martensite and austenite in the SMA 
rod. The frontier of the waves are now more easily identified since abrupt jumps occur in 
the strain distribution, which are caused by phase transformations between austenite and 
martensite. In the displacement evolution now we have only one peak within the simulated 
time span, while two peaks are found in the first experiment. On the three chosen time 
instants, the strain distributions are not as smooth as those at high temperature. At 
t = 0.03 ms, the wave frontier is around x = 0.65 cm, propagating along the positive x 
direction. It indicates that the wave speed is a little lower than that in the first experiment, 
as indicated by the analysis given in section 3. Similarly, there are thermal waves caused 
by the mechanical loading due to the coupling effects. 

For the third experiment, the initial temperature is set at 6' = 210 K. Because only 
martensite is stable at this temperature, the initial condition is chosen such that the SMA 
rod is originally at martensite minus state, for which the displacement is set u = — O.llSx 
so Eq = —0.115. This strain value is one of the local minima of the non-convex potential 
energy plotted in Figure ([T]), and calculated in Equation (fTTj) . Numerical results for this 
case are presented in Figure ([7]). It is seen that the entire SMA rod is divided into two 
domains, one consists of martensite plus (with e ~ 0.115), and the other one - martensite 
minus {e ~ —0.115). The interface between the two domains is driven by the impact 
stress loading, as sketched in the strain evolution plot and the wave profiles on the chosen 
time instants in Figure ([7]). With the given computational conditions, the SMA rod 
is converted from martensite minus to plus from the end which under impact loading, 
the phase boundary is driven to propagate along the negative x direction. Due to the 
hysteretic nature of the phase transformation, the input mechanical energy is dissipated 
continuously during the propagation of the phase boundary, and is not able to convert the 
whole rod into martensite plus. The interface stops at around x = 0.5. Correspondingly, 
there are also thermal waves accompanying the martensite transformation. The wave 
propagation speed is much lower compared to those in the previous experiments, and the 
wave speed changes more remarkably during the propagation process, as indicated by the 
plot of wave profiles plot at chosen time instants. At t = 0.03 ms, the wave frontier is at 
X = 0.5 cm and is unable to move further toward the end a; = 0. 

The forth experiment is to investigate the dissipation effects due to internal friction in 
the material. As analyzed in section 3, the dissipation effects are independent of temper- 
ature, so we set the initial temperature at 6q = 310 K to exclude phase transformation, 
so that the comparison will be easier. The internal friction u is set three times larger 
at 30, and all other computational parameters are chosen the same as those in the first 
experiment. The strain evolution and wave profiles at the same three chosen time instants 
are presented in Figure . By comparing with those in Figure ([H]) , it is observed that the 
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dissipation effects are enhanced, the peak values of the wave decrease faster. At t = 0.01 
ms, the peak value of the wave profile is around 0.075, located around x = 0.6 cm, while 
the counterparts with z/ = 10 are peak value 0.089 at around x = 0.055 cm. This indicates 
that when u is increased, not only the peak value dissipated faster, but the wave speed 
was also slightly slowed down. 

The final experiment is to show numerically the dispersion effects due to the Ginzburg 
term in the wave equation. As indicated by Equation ( fTSj) . the dispersion effects are 
more pronounced at low temperature since ki = ki{6 — 6i) will be smaller. To perform 
the analysis, the initial conditions are set the same with those in the third experiment, 
except that kg is set three times larger at 30. The strain evolution and wave profiles for 
three chosen time instants are presented in Figure By comparing the results with 
those in Figure (JTj), it can be seen that the entire rod is still divided into two domains, 
one for martensite plus {e ~ 0.115) and the other - for martensite minus {e ~ —0.115). 
However, the wave propagation speed is faster with larger kg value. The interface between 
martensite minus and plus is located at around x = 0.35 cm when t = 0.03 ms, while 
for fcg = 10 it is at X = 0.5 cm. This observation agrees with the linearization analysis 
carried out in section 3. At the same time, the entire rod is converted from martensite 
minus to plus, which indicates that smaller amount of input energy is demanded for the 
phase transformation when kg value is larger. In other words, the phase transformation 
becomes easier to deal with when the capillary effects are enhanced. 

From the above numerical experiments follow that nonlinear thermo-mechanical wave 
propagations caused by impact loadings in the SMA rod can be remarkably influenced by 
the material temperature, internal friction, and capillary effects. Thermal waves could be 
induced by impact mechanical loadings. Wave propagation patterns are more complicated 
when phase transformations are involved, and the dynamic response of the material in 
this case is very different from those with no phase transformations. 

6 CONCLUSIONS 

In this paper, a mathematical model for the analysis of wave propagations in a shape 
memory alloy rod induced by a stress impact was constructed. The modified Ginzburg- 
Landau- Devonshire theory was employed for modelling dynamic processes in SMA rods. 
The first order martensite phase transformations and thermo-mechanical coupling were 
incorporated into the model. A multi-domain decomposition method was employed in 
conjunction with the Chebyshev collocation method for spatial discretization, and the 
backward differentiation formula was used for solving the resulting differential-algebraic 
system. The nonlinear thermo-mechanical wave propagations in the SMA rod were sim- 
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ulated with various initial temperatures (with and without phase transformation), the 
effects of phase transformations on the wave propagations were analysed numerically, 
along with the effects of internal friction and capillary. 
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Figure 1: Constitutive relations (left colum) for a shape memory alloy and the associated 
potential energy density (right colum), at various temperatures, (top) 9 — 210 K. (middle) 
e = 240 K. (bottom) 9 = 310 K. 
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Figure 2: Mechanical dissipation due to martensite phase transformation 
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Figure 3: Sketch of domain decomposition and discretization 
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Figure 4: Shape memory alloy rod under impact loadings 
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Figure 5: Nonlinear thermo-mechanical wave propagations in a shape memory alloy rod 
caused by a stress impact, initial temperature is 6* = 310 K. 
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Figure 6: Nonlinear thermo-mechanical wave propagations involving phase transforma- 
tions in a shape memory alloy rod caused by a stress impact, initial temperature is ^ = 240 
K. 
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Figure 7: Nonlinear thermo-mechanical wave propagation in the shape memory alloy rod 
caused by a stress impact, involving martensite phase transformations. Initial temperature 



is e = 210K 
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Figure 8: Numerical analysis of the dissipation effect of internal friction on wave propa- 
gations in a shape memory alloy rod. Initial temperature is6' = 310K, i/ = 30. 
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Figure 9: Numerical analysis of the dispersion effects caused by capillary effects in wave 
propagations in a shape memory alloy rod. Initial temperature is ^ = 210 K, kg — 30. 
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